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Abstract 

A relativistic and quantum mechanical framework to compute nuclear transparencies for pion 
photo- and electroproduction reactions is presented. Final-state interactions for the ejected pions 
and nucleons are implemented in a relativistic eikonal approach. At sufficiently large ejectile 
energies, a relativistic Glauber model can be adopted. At lower energies, the framework possesses 
the flexibility to use relativistic optical potentials. The proposed model can account for the color- 
transparency (CT) phenomenon and short-range correlations (SRC) in the nucleus. Results are 
presented for kinematics corresponding to completed and planned experiments at Jefferson Lab. 
The influence of CT and SRC on the nuclear transparency is studied. Both the SRC and CT 
mechanisms increase the nuclear transparency. The two mechanisms can be clearly separated, 
though, as they exhibit a completely different dependence on the hard scale parameter. The 
nucleon and pion transparencies as computed in the relativistic Glauber approach are compared 
with optical-potential and semi-classical calculations. The similarities in the trends and magnitudes 
of the nuclear transparencies indicate that they are not subject to strong model dependencies. 
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I. INTRODUCTION 



A commonly used variable to map the transition from hadronic to partonic degrees of 
freedom is the nuclear transparency. For a given reaction process, it is defined as the ratio 
of the cross section per target nucleon to the one from a free nucleon. Accordingly, the 
nuclear transparency provides a measure of the attenuation effects of the nuclear medium 
on the hadrons produced in some reaction. A phenomenon finding its roots in QCD is color 
transparency (CT). It predicts the reduction of final-state interactions (FSI) of the pro- 
duced hadron with the surrounding nuclear medium at sufficiently high momentum transfer. 
Thereby, the hadron is created in a point-like configuration (PLC) and propagates as a color 
singlet through the nucleus before evolving to the normal hadron state. If CT effects were 
to appear at a certain energy, the nuclear transparency would be observed to overshoot the 
predictions from traditional nuclear physics expectations. 

Measurements of nuclear transparencies in search of CT have been carried out with the 

A(p, 2p) naaa 

and A(e, e'p) , llOj reactions, p-meson production [Hi, ll2| 

and diffractive dissociation of pions into di-jets [13j. Nuclear transparencies for the pion 
photoproduction process — > ii~p in 4 He have been measured in Hall A at Jefferson 
Laboratory (JLab) 14]. A Hall C experiment has extracted the nuclear transparency for the 
pion electroproduction process ep — > e'n + n in 2 H, 12 C, 27 A1, 63 Cu and 197 Au 15]. Ref. [3] 
reports calculations in a semi-classical model for the latter electroproduction experiment. 
In Ref. 17| . we introduced a relativistic and quantum mechanical model for computing the 
nuclear transparencies for the pion photoproduction reaction and compared its predictions 
to the 4 He(7,p7r~) data and results from a semi-classical model developed by Gao, Holt 



and Pandharipande 



18] . In this paper, we outline the model in more detail and extend 



it to electroproduction reactions. The intranuclear attenuation which affects the ejectiles 
(nucleons and pions) is modeled in terms of a relativistic eikonal approach. The bound-state 
wave functions are obtained from a relativistic mean-field model. At sufficiently small values 
for the de Broglie wavelength, we use a relativistic version of Glauber multiple scattering 
theory. At wavelengths approaching the range of the nucleon-nucleon and pion-nucleon 
interaction length, the model offers the flexibility to use optical potentials for modeling FSI 
mechanisms. Short-range correlations (SRC) induce local fluctuations in the nuclear density. 
These corrections beyond the mean-field approach influence the intranuclear attenuation. 
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The corresponding changes in the nuclear transparencies have been studied in great depth 



within the context of A(e, e'p) reactions 3, 20, 21]. In A(j, Ntt) and A(e, e'Nir) processes, 
both the emerging nucleons and pions are subject to these density fluctuations. The SRC are 
incorporated into our model through the introduction of a well-chosen central correlation 
function which induces density correlations into the final system. In our procedure, the 
proper normalization of the wave functions is guaranteed. 

Section [TT]of the paper presents the formalism used to calculate the nuclear transparencies. 
A factorized expression for the cross section is derived for A (7, Ntt) (III Aj) and A(e, e'Nir) 
( 1IIBI) . Next, in Sect. Ill CI the framework for computing the effects stemming from FSI are 
discussed. Thereby, special attention is paid to a parametrization of the ttN scattering 
parameters which are required in Glauber calculations. In Sect. Ill D I the incorporation of 
the CT phenomenon and SRC is discussed. The results of our numerical calculations are 
presented in Sect. IIHI FSI effects are investigated and transparency results are shown for 
the pion photo- and electroproduction reactions from various target nuclei. Our conclusions 
are stated in Sect. [TV] 

II. FORMALISM 

In this section, the formalism used to describe A(j, Ntt) and A(e, c'Ntt) reactions is 
presented. 

A. Pion photoproduction 

We use the following notations for the four-momenta in the lab frame: q^{q, q) for the 
photon, P^(Ea,Pa = 0) for the target nucleus, P^^Ea-i, Pa-i) for the residual nucleus, 
Pn(EniPn) and P^(E n ,p n ) for the ejected nucleon and pion. The missing momentum p m 
is defined as p m = —pA-i = Pn + Pn — q and the outgoing nucleon has spin m s . The fivefold 
differential cross section in the lab frame reads 

d 5 a MA^mNp^PN f -rsr^\ kA (~f,NTT) 2 m 

free/ j Y^-H 1 V-) 



dE^dVt^dVt N A{2TjfqE A 



with the recoil factor given by 



, Ea-i 

Jrec 



E A 



Ea-i V Pn 
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and A4^' N ^ the invariant matrix element: 

M (y,N*) = p M ms; p^J R M R \6\q^ Pp+) , (3) 

where JrM r are the quantum numbers of the residual nucleus. We restrict ourselves to 
processes with an even-even target nucleus A. 

The wave functions for the bound nucleons are constructed in an Independent Particle 
Model (IPM). We use relativistic wave functions from the Hartree approximation to the 



Walecka-model with the Wl parametrization 221] • For the sake of conciseness, only the 
spatial coordinates of the nucleons are written throughout this work. The single-particle 
wave functions <p a adopt the following form for a spherically symmetric nuclear potential 



231 : 



(j) a (r) = nKm (r, a) 



F nK {r)- 



(4) 



Here, n is the principal quantum number, k and m denote the generalized angular momentum 
quantum numbers. The spin spherical harmonics y± Km are defined as: 

y Km (Q,a) = ^ {lmi-m s \jm)Y lmi (ti)x^ mi> {a) , 



y^ Km {Q,a)= ^2(lm l ^m s \jm)Y Imi ( y Q)xi ms (a) , (5) 



1 K, K> - \ K - 1, K > 

with j = \k\ ,/=< ,/=< 

2 I — K — 1, K < 1 —K, K < 

The ground-state wave function of the target nucleus \P A l + ) = ^/^ s '(fx, . . . ,ta) is obtained 
by fully anti-symmetrizing the product of the individual nucleon wave functions <p a . We 
model the pion photoproduction process by means of a contact interaction: the initial nu- 
cleon, impinging photon, the ejected pion and nucleon, couple in a single space-time vertex. 
As the process can take place on any of the nucleons in the target nucleus, we get the 
following general expression for the corresponding photoproduction operator: 

A 

6 = Y,0,{n). (6) 

i=l 

We assume that O is exempted from medium effects. This is a common assumption in nuclear 
and hadronic physics and is usually referred to as the impulse or quasi-free approximation 



(IA). In the context of A (e, e'p) reaction, for example, the impulse approximation provides 
a fair description of the data [24\. It is also applied in the experimental analysis of Ref. [ijj] 



and the model of Ref. [25( . The impinging photon with polarization A is represented by 

A^(X, = e^(A)e^ . (7) 

Here, e M (A) is the polarization four- vector of the photon. The wave function of the ejected 
nucleon is written as 

\ P N m s) = *^ mj ifi) = & N , N {?i\ n,..., f j¥i , . . . , f A )u(p N} m a )e#»* , (8) 

which is the product of a positive-energy Dirac plane wave <j>p N and an operator S N , N . This 
operator describes the attenuation of the ejected nucleon through soft final-state interactions 
with the other nucleons. The wave function for the ejected pion adopts a similar form as 
the nucleon one, i.e. a plane wave convoluted with a FSI factor S\ N : 

= 4t } ft) = • • • > • • • > ■ ( 9 ) 

The final A-nucleon wave function is constructed by anti-symmetrizing ^j^ m with the wave 



function for the residual nucleus i $ J £^ a \ 



A 



^{rl;^, . . ,f A )«(p iV ,m s )e*- fl #^, ■ • • • (10) 



As \l/^ s ' and \fr p ^' ms are fully anti-symmetric, each term of the operator ([6]) will yield the 
same contribution to the matrix element ([3]) and we can restrict ourselves to the term with 
coordinate r*i and multiply it with A. With the above expressions for the operator and the 
wave functions of the hadrons involved in the reaction, we can write for the matrix element 
of Eq. (EJ) in coordinate space: 



M (j,Nir) =A I ^ I I tf A 



t 



x e-^&Jv(ra;F 2 ,...,fAW • (H) 

We assume that S^'n and S^n are spin independent and that only elastic and mildly inelastic 
collisions with the spectator nucleons occur. The actual nuclear transparency measurements 
select events whereby the undetected final state with (A — 1) nucleons \P^_ v JrM r ^ is left 



with little excitation energy, which makes these assumptions very plausible. In computing 
the matrix element of Eq. (fTTj) we consider processes of the type displayed in Fig. [IJ The 
following spectator approximation is assumed to be valid for a struck nucleon with quantum 
numbers ol\ : 

t 



dfi ... / dr 



0^(P n (r 1 ))«S^ JV (P ri (f 1 ); P n (f 2 ), . . . , P„(f A ))0 Q2 (P n (f 2 )) . . . <f> aA (P n (f A )) 

xe-^' rl 4iv(n; r 2 , . . . , r j4 )O M (f 1 )e i ^ 1 ai (P m (f 1 ))0 Q2 (P m (f 2 )) . . . QA (P m (f A )) 

~ h n {r 2 )p m {r 2 ) ■ ■ ■ h n {r A )p m {? A ) J dn . . . J df^t^ (f^S^'Tv (ri; P„(r 2 ), . . . , P„(r A )) 

e-^^^Cfi; r 2 , . . . ,fA)O^(r 1 )e if ^0 ai (P m (r 1 ))|0 a2 (P m (r 2 ))| 2 . . . |0 aA (P m (r4))| 2 , 

(12) 

with P m and P n permutations of the set {f*i, . . . , r^} occurring in the anti-symmetrization of 
the nucleon wave functions. Due to the presence of the delta functions, the rhs of Eq. f[T2"j) 
is non-vanishing under the condition that P m {fi) = fi and P m {fi) = P n (r*j) for z = 2, .., A 
This means that both the bound wave function ot\ and the ejected nucleon have the same 
spatial coordinate as the operator, f\. Moreover, all (A — 1)! permutations of the subset 
{r 2 , . . . , f A } yield an identical rhs. 

Thus, after expanding the wave functions in Eq. ( fill and employing Eq. ( fl2l) . we arrive 

at 



ifo.Nir) ^4(^4 — 1) 

[}■ ~ J ar x J ar 2 ... J ar A 



'0.1 



{f 2 )\\ . .\<j> aA {f A )f 



(13) 



xu\p N , m s )S nN (r 1 ; r 2 , . . . , rx)<SjV'jv(fi; r 2 , . . . , rA)e M (A)O M (r*i)e" iPm ' ri ai (fi) 

We now define the FSI factor J-fsi(^) : 

-pFsi(r) = y dfi... J df A \(p a2 (f 2 )\ 2 . . . \(p aA (fA}\ 2 <SnN(f;f2, ■ ■ ■ ,fA}S N >N(f; r 2) ...,f A ), 

(14) 

and write 

M% Nn) « / ^FSi(ri)« t (piv,m s )e^A)0 M (f 1 ) e - i ^-V ai (r 1 ). (15) 

In what follows, we assume that the pion production operator acts on a bound-state wave 
function as a scalar (factorization assumption): O M (r)0 ai (r) = C(j) ai {f). With 

< (P) = y tfV-*^ (rl-P FSI (r) , (16) 
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a, 



> ^^iTiiiil 

























a 



A 



FIG. 1: (Color online) Diagram included in computing the matrix element of Eq. (Ilip . The dashed 
lines denote the FSI of the ejected pion (red) and nucleon (blue) with the spectator residual 
nucleons. The diagram shown here is representative for the spectator approximation: one active 
nucleon N and tt are subject to soft collisions with frozen spectator nucleons which occupy the 
single-particle levels ct2, 03, . . . , cxa and are not subject to changes in their quantum numbers. 



we can write 



(17) 



When studying nuclear transparencies, it is convenient to factorize the invariant matrix 
element such that it becomes a convolution of a factor describing the elementary pion pho- 
toproduction process and a factor modeling the combined effect of all FSI mechanisms of 
the outgoing hadrons. To reach this goal we relate the 7 + A — > (A — 1) + N + ir matrix 
element in Eq. ( TPTj) to the one for for free nucleons 7 + iVj — > N + 7r 



u ] (p N , m s )e fl (\)O IJ ,(f 1 )u(p m , m s >) 



(18) 



with m s i the spin of the initial nucleon. First, we consider the situation with vanishing 
FSI, second the more realistic case with inclusion of a FSI phase operator. When ignoring 
FSI, the wave functions for the ejected hadrons reduce to plane waves and ^fsi(^*) = 1> 
4>ai(Pm) = 0Qi(Pm)- After substituting in Eq. (jTTI) the completeness relation for Dirac 



7 



spinors: 

[u(p m , m' s )u(p m , m' s ) - v(p m , m' s )v(p m , m' s )\ = l 4x4 , (19) 

one obtains 

V J / RPWIA ^ — ' V J / m s ,m' s 

m' s 

— negative energy terms , (20) 

where the RPWIA denotes the relativistic plane wave impulse approximation. From this 
last expression it is clear that even with vanishing FSI the presence of negative-energy 
components makes factorization impossible. In what follows we neglect those terms: 

(M^LwT. * (M%£?) , uiPmM^ipm). (21) 

V J /RPWIA A — ' V / m s ,m' s 

m' s 

The contraction of the Dirac spinor u with the bound nucleon wave function <p ai if negative 
energy components are neglected is given by 



u(p m ,m' s )(f) ai (p m ) = — ^ J" — "^(PmlXi^J^ffip, a) , (22) 



where m N . is the free mass of the bound nucleon, E N .{p m ) = \Jrn 2 N . + pj n and 

Ow(Pm) = T-i 2mNl 9nn(Pm)- (23) 

h, N . + m N . 

In this last equation g nK is defined as 

9nM =iJ^T r 2 dr^l 3l (pr) , (24) 
V vr J r 

with ji(pr) the spherical Bessel function of the first kind. After squaring the matrix element 
and summing over the quantum number m of the bound nucleon wave function, one can use 
the following property of the spin spherical harmonics y Km 

Y.y^ P ,d)yi m {n p ,o) = (^±±)i 2x2 . (25) 

m 

Finally, by using x\ X±m> = &m s m' , the free pion production process can be formally 

2 i m s 2 ' s s 

decoupled from the typical nuclear effects: 

}i A,m,m s 



X 4 

X,m s ,m 
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The matrix element for the pion production process on a free nucleon is related to the cross 
section by 

- V I [m^ N7T A \ 2 = 4?r ( g - m U 2 d ^ n (27) 

4 I V / ifree J ms ,m' s 1 m N .m N d\ty K ' 

\,m a ,m' s 1 

with s = {p^ N + p%) 2 and t = (g M — p%) 2 the Mandelstam variables of the free process. 

After substituting Eqs. (1261) and (1271) in Eq. ([1]), the differential cross section for 7 + A — > 
(A — 1) + iV + 7r in the relativistic plane wave impulse approximation (RPWIA) reads 

d 5 a \ M A _ l p n p N (s - m 2 J 



dE^dVt^dVtN J RPWIA 4irm Ni qE A 



' J 1 



rcc 



x 2 3 + l(E Nt (p m ) + m Ni ) 

4tt 2m N . 1 nKKymJl d\t\ K ' 

When FSI are included, the derivation outlined earlier is no longer possible due to the 
presence of .Ffsi(^) i n 4>a- We define a distorted momentum distribution along the lines of 



Ref. [26| 



pD(Pm) = ^ l M (^' m s)^fi(Pm)| 2 - ( 29 ) 

m a ,m 

When FSI and negative energy contributions to 0„ are neglected, Eq. (129]) reduces to 
2mx. mNi 1 a nn (Pm) | 2 - Based on this analogy, we write the differential cross section 
with FSI as 

( d 5 a \ M A _ lPnPN (s-m 2 N f ^ _ da^ 

\dE n dn n dQ N ) D ^ 4nm Ni qE A JrecPniPN^mjjjjj ■ {W) 

B. Pion electroproduction 

The four-momentum of the virtual photon 7* is q^{u,q) and the z-axis lies along q. The 
incoming (scattered) electron has four-momentum p^(E e ,p e ) (p^,(E e i,p e >)) and spin s (s'), 
9 e denotes the electron scattering angle. With these additional notations and conventions, 
the differential cross section in the lab frame reads 

(Pa m\p e , M A _ x m N p^p N I A ^(7*,iV7r) 2 , qi v 

free \ M fi ' V^ 1 / 



dVteidE^dE^dVt^dVLN (2vr) 3 p e 2(2nfE A 

with the recoil factor f rec as in Eq. (j2J). The invariant matrix element Ad^ eN can be 
written as 

M (e,e>N,) = (p , P ^ ljR M R \j^\ P p+) , (32) 



with the electron current 



3n = U(p e >,s')-f^,u(p e ,s) 



Q 2 = —q^q^ and the hadron current J**. By defining an auxiliary current 

Jo 

S = j m 

and using current conservation, the following identity can readily be proved: 

= -ciiJi = -didijJj = - a i e i (A)e i (A) Jj , 

x={x,y,z) 



(33) 



(34) 



(35) 



where e(A) is the unit vector along the axis A = (x, y, z). After defining the electron density 
matrix 

Pay = £ ■ $ WW ■ «] ( 36 ) 

ss' 

and the hadronic matrix elements 



w x = (P^P^P^JrMrI^X) ■ J\P£0 + ) 



we can write for the matrix element 

EK 

ss' 

With the degree of transverse polarization defined as 



( {e,e'Nir) 
fi 



W XX' 



2q 2 2 6 e 



-i 



27] 



the electron density matrix becomes 

( 1(1 

Q 2 l 

PXX' 



mi 1 — e 







V -| A /2^e(l + 6) 







2ge(l + e) 


02 



\ 



/ 



(37) 



(38) 



(39) 



(40) 



After substituting Eq. (jlQl in Eq. (!38|) . one can factor out a part containing all the variables 
related to the electrons in the differential cross section: 

d s a _ ci 5 cr„ 



d£l e >dE e idE n d£l n dQ,N dE^dVL^dVL 



N 



iV7r|2 



(41) 



hhrc /vlir^ = (P^,P;m s ,P^_ 1 J i? M R |a|g^,P^0 + ), C = %fflgg" and r 



EjRy 1 



2^--^-^— is the electron flux factor, with the virtual photon equivalent energy P 7 = -^mj 
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the fine-structure constant a, and s = (g M + P^) 2 the Mandelstam variable of the virtual 
photoproduction process. The cross section can be cast in the following form 

d 5 a v d 5 a T d 5 a L d 5 a TT rr—^ d 5(T TL 

~ + e , „ — — he \-y/e(e+l) 



dE^dVL^dVL^ dE^dVL^dVL^; dE^dVt^dVtis! dE^dVL^dVL^ dE^dVt^dVt^ 

(42) 



with 



S " T c E IM+W 



dE^dVt^dVt^ 2 

wis -Mr 



d 5 a L O 





Z I 5 



^ E 



dE^dVL^dVtM ijj 

m s M R 



<r ' air c £ 0^i 2 -i^i 2 ] > 



dE^dVt^dVtN 2 

m s M H 



im---4fS^W,,r M . (43) 



As for the photoproduction case, we wish to establish a relation between the invariant 
matrix element for virtual-photon pion-production on a nucleus (AiJ^^) and on a free 
nucleon {■M'j i , ^ e ). In comparison with the real photoproduction process, the virtual photon 
has an extra degree of polarization and Q 2 ^ 0. This does not alter the derivation presented 
in the previous subsection and after neglecting negative energy contributions, one arrives at 

Mf^ « (2nf 2 J2( M }lOx,rn s ,rn 3 MPrn, M<^(pJ ■ (44) 

The matrix element A4X' fre g is related to the free electroproduction process by 

d 5 a eN _ r , mj 

dE e ,dQ e ,d0*d | t | ~ 2(27r) 2 (s' - m^, 



r '^fc^Ei^?:£i 2 . («) 



where V = -Tr^^-^-r— is the electron flux factor, with the virtual photon equivalent energy 

/ 2 

K = s ~ m y . Further, s' = (p^ + p^) 2 and i = (g M — p^) 2 are the Mandelstam variables for 
the free process. Starred variables denote center-of-mass values. 

With pd defined in Eq. ( T29l) and by making use of Eqs. (H4l) and (H5]L we arrive at the 
factorized form for the differential A(e, e'Nir) cross section: 

rfV ^ rM,_ M ( S '-m^) 2 x ci 5 a eiV 

-JrecPD , 7177711 ' l 4fc V 



dQe'dEe/dE^dQ^dQN y D P 2m^E^EA rec dE e >dQ e id \ t 

We wish to stress that the assumptions made to arrive at this expression, are essentially 
identical to those made for the real photon case discussed in the previous subsection. 
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C. Final-state Interactions 



The Glauber approach can be justified when the wavelength of the outgoing hadron is 
sufficiently small in comparison to the typical interaction length with the residual nucleons. 
In the context of A(e, e'p) reactions 28] it was shown that the Glauber model represents 
a realistic approach to FSI for proton kinetic energies down to about 300 MeV. This cor- 
responds to proton de Broglie wavelengths of the order of 1.5 fm. For pions comparable 
wavelengths are reached for kinetic energies of the order of 700 MeV. 

ativistic Multiple-Scattering 



24]]. In the RMSGA, the wave 



A relativistic extension of the Glauber model, dubbed the Re 
Glauber Approximation (RMSGA), was introduced in Ref. 
function for the ejected nucleon and pion is a convolution of a relativistic plane wave and an 
Glauber eikonal phase operator which accounts for FSI mechanisms. In Glauber theory the 
assumption is made that a fast moving particle interacts through elastic or mildly inelastic 
collisions with frozen point scatterers in a target. Scattering angles are assumed small and 
each of the point scatterers adds a phase to the wave function, resulting in the following 
expression for the Glauber eikonal phase: 

A 

S m (r, r- 2 ,...,r 4 ) = n[l" 1^(6- bj)6( Zj - z)] (with i = vr or N') . (47) 

i=2 

Here, r*j(bj, Zj) are the coordinates of the residual nucleons and r(b, z) specifies the interaction 
point with the (virtual) photon. In Eq. (l4"Tj) . the z axis lies along the path of the ejected 
particle i (the proton or pion), b is perpendicular to this path. The Heaviside step function 
9 guarantees that only nucleons in the forward path of the outgoing particle contribute to 
the eikonal phase. 

Reflecting the diffractive nature of the nucleon-nucleon (N'N) and pion-nucleon (wN) 
collisions at intermediate energies, the profile functions T N > N and T nN in Eq. (T47j) are 
parametrized as 

T iN (b) = exp (~) (with i = n or N>) . (48) 

Here, the parameters (total cross section), [3iN (slope parameter) and (ratio of the 
real to imaginary part of the scattering amplitude) depend on the momentum of the outgoing 
nucleon or pion i. For i = N', we determined the parameters by performing a fit 24j to the 
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N'N — > N'N databases from the Particle Data Group (PDG) |29|. For the pion, q£ 



fitted to data collected by PDG 



29j . The analysis of the slope parameter in Ref . 30] was 



r tot 
riV 



was 



used for the f3 n N fits. Fits provided by SAID 
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32] and data from PDG 29J were used in 
and of Figs. [2131 and H] are the result 



constructing the fits for e^jv- The fits for cr^,/^ 
of a x 2 minimization of the data against a a n-th degree polynomial (with n < 10). An 
alternative way of determining f3 n N, is via the relation 



01 



(J 



tot \2(^ _|_ £ 2 



ttN , 



ttNj 



N 



(49) 



with cr^y the elastic cross section. Fits for cr^y to data from PDG are also presented 
in Fig. [2j The two sets for the (3^ parameter in Fig. [3] do not produce significantly 
different results for the numerical calculations presented here. We use the % 2 fit for fi nN in 
all calculations presented in this paper. 

The Glauber operator of Eq. (PIT]) is an A-body operator. As a consequence, it requires 
integrations over all spectator nucleon coordinates in Eq. ffT4l) . which is computationally 
very demanding, in particular for heavy target nuclei. In 7W + A — > (A — I) + N + n 
calculations, a product of two Glauber phases is involved and the cylindrical symmetry of 
the individual phases is lost. A Romberg algorithm is used to perform the integrations over 
the spatial coordinates in Eq. (1T41) . 

For nucleons with a kinetic energy lower than about 300 MeV, the approximations under- 
lying the Glauber formalism are no longer applicable, and an alternative method to model 
FSI is required. Under those circumstances our framework provides the flexibility to adopt 
the Relativistic Optical Model Eikonal Approximation (ROMEA) [33|. In the ROMEA ap- 
proach, the wave function of a nucleon with energy E - 



+ m? N after scattering in a 
scalar {V s (r)) and vector {V v (r)) spherical potential has the following form: 



(f) 



E + m n 



2m 



N 



E+m N +V s (r)-V v (r) 



a ■ p 



(50) 



with the eikonal phase determined by 



iS N > N (b,z) 



*1T 



dz 



V c {b, z') + V so {b, z') a ■ (b x K) — iKz' 



(51) 



In this last equation, K = |(fcj + kf) is the average of the initial and final momentum of 
the scattering particle. In the small angle approximation, K and points along the 
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p [GeV/c] 10 




pJGeV/c] 10 



FIG. 2: The pion lab-momentum dependence of the data 29| and adopted fits for the total and 
elastic cross section for ir~ — p (upper panel) and tt + — p (lower panel) scattering. 

z-axis. The central and spin-orbit potentials V c and V so are functions of V s and V v and their 
derivatives 



Additional approximations were used in the implementation of optical-potential FSI in 
this ROMEA model. The dynamical enhancement of the lower components of the scattering 
wave function (1501) is ignored as at low momenta the lower components are small compared 
to the upper components due to p and at higher momenta (V s — V v ) is small in comparison 
to (E + mjy). The operator p was also substituted by the asymptotic value pW- Finally, as 
collisions were assumed spin- independent in ( 1T21) . the spin-orbit potential V so in Eq. ( 1511) is 
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p [GeV/c] 10 




p [GeV/c] 10 



FIG. 3: The pion lab- momentum dependence of the data 30(| and fits for the (3p n parameter for 



7r — p (upper panel) and 7r + — p (lower panel) scattering. Full curves are a x 2 fit to the data, 
whereas the dashed curves result from Eq. (|49p . 



(52) 



neglected. This yields the following phase factor entering in Eq. ([HI) : 

S N > N (r) 



In contrast to the Glauber eikonal phase, the optical potential eikonal phase of Eq. (I52p 
depends solely on the coordinate r which defines the interaction point. As a consequence, 
it can be taken out of all the integrations in Eq. (1141) and the cylindrical symmetry of the 
pion Glauber eikonal factor is retained, hereby considerably reducing the cost of computing 
the total FSI factor ^Vsi- For the numerical evaluation of the ROMEA phase factor, we 
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p [GeV/c] 10 




p [GeV/c] 10 



FIG. 4: The pion lab-momentum dependence of the ratio of the real to imaginary part of the ir p 
(upper panel) and ir + p (lower panel) amplitudes. The diamonds represent an analysis of the data 



by the Georges Washington University group 
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32j, whilst the solid circles are from PDG 



20]. 



The solid line is the fit to the data which is used in the numerical calculations. 



made use of the optical potentia 
parametrization of Cooper et al. 



35] 



van Oers et al. 34( for 4 He and the global (S — V) 
for heavier nuclei. 



D. Color transparency and short-range correlations 

We implement color transparency effects in the usual fashion by replacing the total cross 
sections <r-^ in the profile functions of Eq. (l4"8j) with effective ones [361 ] . The latter induce 
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some reduced pion-nucleon and nucleon-nucleon interaction over a typical length scale lh 
corresponding with the hadron formation length (i = it or N') 



tot 
a iN 



Z < n 2 k 2 > ( Z 



eff 

h H V 1 " h ) \ 9{lh ~ Z) + 9{Z ~ k) ) • (53) 

Here, n is the number of elementary fields (2 for the pion, 3 for the nucleon), k t = 
0.350 GeV/c is the average transverse momentum of a quark inside a hadron, Z is the 
distance from the interaction point and lh — 2p/AM 2 is the hadronic expansion length, 
with p the momentum of the final hadron and AM 2 the mass squared difference between 
the intermediate prehadron and the final hadron state. We adopt the values AM 2 = 1 GeV 2 
for the proton and AM 2 = 0.7 GeV 2 for the pion. TC is the hard scale parameter that governs 
the CT effect. It equals the momentum transfer t = (g M — p%) 2 (pion CT) or u = (q^ — p^ N ) 2 
(nucleon CT) for pion photoproduction and Q 2 for pion electroproduction. Fig. illustrates 
the predicted difference of the CT effect on the pion-nucleon and nucleon-nucleon effective 
interaction. Reflecting its mesonic nature, the pion has a longer formation length and during 
its formation its interaction cross section with the residual nucleons is more strongly reduced 
than for a nucleon. 



0.5 



Z[fin] 



FIG. 5: Comparison of the CT effect on the total effective cross section a^fj for nucleon-nucleon 
(full) and pion-nucleon (dashed) interactions. We consider the situation whereby the ejectile pos- 
sesses a lab-momentum of 2.5 GeV/c. For the hard-scale parameter we adopt 7i = 1.8 (GeV/c) 2 . 



We now proceed with introducing a method which allows us to implement the effect of 
SRC in the relativistic Glauber calculations. The proposed method adopts the thickness 
approximation as a starting point. In the thickness approximation, the density |p ai (r*j)| 2 of 
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the individual nucleons in Eq. (JTlj) is replaced by an averaged density (r) defined as 

p [l A \f) = A J df 2 ... J dr A (^\f,f 2 ,...,r A )) i ^-(r,f 2 ,...,r A ). (54) 



In terms of p A (f) the FSI factor of Eq. (fT4"l) can be approximated by 
•^FSi Ck ( f ) = -JaZi I df 2 ... I df A p [ ^ ] (f 2 )p [ ^ ] (f 3 ) . . . pjj ] (f A )S nN {f; f 2 , . . . , r A )S N > N {r; r 2 , . . . , f A ) 



(55) 



In combination with the operators of Eq. (H7|) the expression can be further simplified to 



j^it (f) 



dr- 



'' A 



1 — ^N'p(bN' — bN'2)Q{zN'2 — Zn>) 



x dr 



1 — F n p(b n — b n2 )9(z 7r2 — z w 

p [l M) 



A 



1 — ^N'nipN' — bN'3)8[ZN'3 ~ Zn>) 

1 — ^wn(b n — ^3)^(^3 — z n) 



(56) 



where t z is the isospin (1 for protons and -1 for neutrons) of the nucleon on which the 
initial absorption took place. The zn> (z w ) axis lies along the ejected nucleon (pion). The 
above expression is derived within the context of the IPM. It is clear that the nucleus has 
a fluid nature and that the IPM can only be considered as a first-order approximation. In 
computing the FSI effects by means of the Eq. (|56|) one fails to give proper attention to 
one important piece of information: namely that one considers the density distribution of 
nucleons given that there is one present at the photo-interaction point r. 

The two-body density p A (ri, f 2 ) is related to the probability to find a nucleon at position 
r*2 given that there is one at a position f\. We adopt the following normalization convention 
for Pa 

J dn J df 2 p [2 A ] (ft, r 2 ) = A(A - 1) . (57) 

In the IPM on has (p A (fi,f 2 )) = ^^P^ A {fx)p^ A (f 2 ). The nucleus has a granular struc- 
V / IPM 

ture as the nucleons have a finite size. This gives rise to strong nucleon-nucleon repulsions 
at short internucleon distances which reflect themselves in short-range correlations (SRC) 
at the nuclear scale. One can correct [p A (fi, f 2 ) J for the presence of the SRC by adopting 
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the following functional form 37 1 



pA ] (n,f 2 ) = 7(fi) (p^in,^)) 7(^2)^(^12) = — - ^l{ri)p [ A{ri)p { A{r2)l{r2)g{ri2) , 

V / IPM A 

(58) 

with g(r 12 ) the so-called Jastrow correlation function and 7(f) a function which imposes the 
normalization condition of Eq. (jSTj) on pj^ (ri,r 2 ). The function 7(f) is a solution to the 
following integral equation 



7(r*i) / rfr 2 pJJ 1 (r 2 )5f(r 1 2)7(f 2 ) = A. 



(59) 



which can be solved numerically. The Glauber phase factor of Eq. ( 156]) can now be corrected 
for SRC through the following substitution 

A p [ a (r 2 ,r) 



P2 ] (r 2 ) 



^4-1 pW(r) 



7(^2)^2 (^2)7(0^(1^2 - r]) = pf(r 2} r) 



(60) 



whereby pY(r2,r) adopts the expression fl58l) . These manipulations amount to the following 
final expression for the Glauber FSI factor including SRC: 



7(r)2p2 1 (r 2 )7(r)^(|r2 - r]) 



.4 



— * — * 



Tz+l 



1 — Fnp^bir — b^2)Q[z^2 ~ Z n ) 



X 



7(^2Pa if'- 2)1(^)9(1^ i - A) 



dr- 



A 



1 — ^N'n{bN' — t>N'?>)9{zN'-Z ~ Zn>) 



1 — Tvi-n^Tr — b^)6{z^ — z n ) 



■ (61) 



The effective density of Eq. (1601) accounts for the fact that the motion of each nucleon does 
depend on the presence of the other ones. In Fig. [6] we display the effective nuclear density 
as it would be observed by a nucleon or a pion created after photoabsorption on a nucleon 
at the center of the nucleus. The figure shows the density for Fe as computed in the IPM 
(p^(x,y, z = 0)) and with the expression based on the substitution of Eq. ( 1601) 

lix,y,z = 0)p [ A ] (x } y,z = 0)7(2; = 0,y = 0,z = 0)g(\r\) . 

In Fig. p and all forthcoming numerical calculations we use a correlation function g(\r\) from 



Ref. 



38| . It is characterized by a (Gaussian) hard core of about 0.8 fm and a second bump 
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which extends to internucleon distances r of about 2 fm and reaches its maximum for r\% 
1.3 fm. This correlation function provided a fair description of the SRC contributions to 
12 C(e, e'pp) 39] and 16 0(e, e'pp) 



40 1 . It is clear that the SRC lead to a local reduction - with 



size of the nucleon radius - of the density around the nucleon struck by the (virtual) photon. 
In order to preserve the proper normalization, this reduction amounts to some enhanced 
density at distances of about twice the nucleon radius. With regard to the intranuclear 
attenuation, the reduction of the density in the proximity of the struck nucleon will result in 
some enhanced transparency close to the photo-interaction point r. The enhanced density 
at positions of about twice the nucleon radius from the struck nucleon, can be expected to 
have the opposite effect. 



III. NUMERICAL RESULTS 



A. The FSI factor 



In this subsection we present a selected number of results of the numerical calculations 
of the RMSGA FSI factor of Eq. (TH1) . We consider the 12 C(7,£>7r~) reaction in a reference 
frame with the z axis along the momentum p^ of the ejected nucleon and the y axis along 
Pn x Pn- The coordinate f denotes the interaction point with the external photon. The FSI 
factor is plotted versus the spherical coordinates in this frame. 

In Fig. [7J we present the calculated norm and phase of the FSI factor in the scattering 
plane (0 = 0) for p^ ~ 2.6 GeV and p^ ~ 2.3 GeV, which are conditions for which Jefferson 
Lab collected data. We present the FSI factor for the proton and the pion separately as well 
as the combined effect when the two are detected in coincidence. 

When looking at the 9 dependence, it becomes clear from Fig. [7J that the norm is smallest 
in the direction opposite the momentum of the particle (being 180° for the nucleon and 
180° — O^tt for the pion). For these directions and large r, the nucleon or pion is created 
close to the surface of the nucleus on the opposite side of its asymptotic direction and has 
to travel through a thick layer of nuclear medium before it reaches a free status. As for the 
r dependence, we see for the nucleon a reduction of the FSI effects for rising r at angles in 
the neighborhood of 9 = 0°, respectively an increment for rising r at 9 = 180°. This is again 
due to the fact that the outgoing nucleon traverses less, respectively more nuclear matter on 
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FIG. 6: (Color online) The effective nuclear density p e ^{r2,r) at Z2 = for He (left) and Fe (right) 
before (upper) and after (lower panel) the inclusion of SRC effects. The effective nuclear densities 
here refer to the situation whereby the (virtual) photon is absorbed at the origin (x = 0, y = 0, z = 
0). 

its way out of the nucleus. The same observations apply for the pion, albeit at the angles 
6Vtt an d 180° — 9^ n . The total FSI factor combines the intranuclear attenuation effects on 
the nucleon and pion. Hence, the norm shows the largest reduction at 9 around 180° and 
180° — 9ntt- The phase of the FSI factor exhibits similar behavior, with the largest phase 
shifts occurring at the discussed angles. 

Fig. [8] teaches us a couple things about the (ft dependence of the FSI factor. As the 
outgoing nucleon lies along the z axis there is no dependence on the azimuthal angle because 
of the cylindrical symmetry. Again, we can see that the absorption is largest when large 
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FIG. 7: (Color online) Radial and polar-angle dependence of the norm (left) and phase (right) 
of the FSI factor .Ffsi i n the scattering plane (4> = 0°) for the 12 C(7,p7r _ ) reaction from the 
lsi/2 level. For the upper (middle) panels, solely the FSI effects on the ejected proton (pion) are 
considered. The lower panels include the net effect of both the pion and nucleon FSI effect. The 
results are obtained for p N = 2638 MeV, p w = 2291 MeV, 6 Nw = -65.19°. 
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amounts of nuclear matter need to be traversed (i.e. large 9). Looking at the pion we see the 
largest attenuation occurs in the upper hemisphere (cos0 > 0) as a pion that is created in 
this region has to traverse the inner core of the nucleus. The combined effect of the pion and 
nucleon contributions is contained in the bottom panel. As the reaction takes place in the xz 
plane, the total FSI factor retains the following symmetry: J-Fsi( r , 4>) — •^ r Fsi( r , @, 27r — 0). 



B. Pion photoproduction 

The experiment E94-104 at Jefferson Lab extracted nuclear transparencies for 7+ 4 He — > 
p + ir~ + 3 He. The measurements were performed for photon energies 1.6 < q < 4.2 GeV 
and for center-of-mass angles 9 c ni = 70° and 90°. In total, the nuclear transparencies 
were measured for eight kinematical settings. In a proposal for a follow-up experiment, 
seven additional kinematics are suggested for measurements at higher photon energies and 



$c.m. = 90° 4lj. We have performed calculatios for the completed and planned experiments. 
Table [I] provides a list of the kinematics. 

We aim at performing calculations which match the kinematic conditions of the experi- 
ment as closely as possible. We use the following definition for the transparency: 

Y. a ! dqY(q) $ dp r , 



rp _ V"^^" - ""^"""^ / RMSGA ^g2) 

£ Q / dqY{q) J dp m ( dE Jn: id n N \ pwiA 
The integrations J dq J dp m in Eq. (1621) were evaluated with a random integration algorithm. 
To this end, random events within the photon beam energy range, detector acceptances and 
applied cuts for each data point were generated for the calculation of the transparency until 
convergence of the order of 5% was reached. Typically, this involves about a thousand events 
for each data point. In Eq. ( 1621) . J2 a extends over all occupied single-particle states in the 
target nucleus. All cross sections are computed in the lab frame. Y(q) provides the weight 
factor for the generated events. It includes the yield of the reconstructed experimental 
photon beam spectrum [ijj] for the photon energy of the generated event. We assume that 
the elementary 7 + n — > n~ + p cross section in Eqs. f )28p and (1301) remains constant 
over the kinematical ranges j dq J dp m which define a particular data point. With this 
assumption the cross section ^2Z1 cancels out of the ratio (16"2"I) . For all kinematic conditions 
of Table [U the pion and nucleon momenta are sufficiently high for the RMSGA method to 
be a valid approach for describing the FSI mechanism. 
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FIG. 8: (Color online) Polar- and azimuthal- angle dependence of the norm of the FSI factor J-fsi 
at a distance r = 3 fm from the center of the nucleus for the 12 C(7,p7r~) reaction from the lsi/2 
level. Separate contributions from the nucleon (upper panel) and the pion (middle panel), as well 
as their combined effect (bottom panel) are shown. Kinematics as in Fig. [7J 
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"cm. 
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Pn 


t) n 


1648 


70° 


989 


47.39° 


1238 


-36.02° 


1648 


90° 


1277 


37.37° 


1015 


-47.73° 


2486 


70° 


1322 


44.37° 


1794 


-31.02° 


2486 


90° 


1740 


34.45° 


1438 


-43.18° 


3324 


70° 


1642 


41.74° 


2363 


-27.56° 


3324 


90° 


2195 


32.01° 


1866 


-38.57° 


4157 


70° 


1949 


39.51° 


2929 


-25.05° 


4157 


90° 


2638 


30.01° 


2291 


-35.18° 


4327 


70° 


2011 


39.1° 


3044 


-24.6° 


4327 


90° 


2727 


29.6° 


2377 


-34.6° 


5160 


70° 


2307 


37.3° 


3606 


-22.8° 


5160 


90° 


3161 


28.0° 


2797 


-32.1° 


6059 


70° 


2622 


35.6° 


4211 


-21.2° 


6059 


90° 


3625 


26.6° 


3250 


-29.9° 


7025 


70° 


2956 


33.9° 


4861 


-19.8° 


7025 


90° 


4120 


25.2° 


3735 


-28.0° 


8057 


70° 


3309 


32.4° 


5555 


-18.6° 


8057 


90° 


4646 


24.0° 


4253 


-26.3° 


9156 


70° 


3683 


31.0° 


6294 


-17.6° 


9156 


90° 


5204 


22.8° 


4805 


-24.8° 


10322 


70° 


4077 


29.7° 


7077 


-16.6° 


10322 


90° 


5794 


21.8° 


5389 


-23.5° 



TABLE I: Central values for the photon energy (MeV), proton momentum pj^ (MeV), proton angle 
#jv, pion momentum p w (MeV) and pion angle 9 n for 6* c m . = 70°, 90°. Angles are measured relative 
to the incoming photon momentum. 
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For a discussion of the computed results compared to the experimental data and a semi- 



classical model we refer the reader to Ref. 



171 ] . In Fig. [9] the separated transparencies for 



the outgoing proton and pion are displayed next to the full result. It is clear from this figure 
that the rise of the transparency at low \t\ can be attributed to the proton contribution. 
This rise can be attributed to the local minimum in the total nucleon-nucleon cross section 
for nucleon momenta of about 1 GeV 

Fig. [9] also shows that the 4 He nucleus is more transparent for pion emission than for 
proton emission. This can be partially attributed to the lower pion total cross sections. As 
pointed out in Fig. [5] the larger formation length, and corresponding bigger reduction of the 
effective cross section make that the CT effect is larger for pions than for protons. In Fig. [TU] 
the computed increase in the nuclear transparency caused by CT and SRC mechanisms 
is shown as a function of \t\. One observes that SRC mechanisms increase the nuclear 
transparency by about 5%. As there is no direct dependence on the hard scale, the increase 
is almost independent of \t\. The CT phenomenon, on the other hand, shows a linear rise 
from almost to over 20% at the largest values of \t\. For —t < 2.5 GeV 2 the predicted 
effect of SRC is larger than the increase induced by the CT mechanism. The SRC decrease 
the slope in the — t-dependence of the CT phenomenon. Indeed, the SRC induces holes in 
the nuclear density in the direct neighborhood of the interaction point (see Fig. [6]) where 
the CT effects are largest. At high \t\ the short-range correlations have a modest impact on 
the magnitude of the CT effects. Our investigations show that by studying the hard scale 
dependence of the transparency the CT-related mechanisms can be clearly separated from 
the SRC ones. 

In the search of phenomena like CT in transparency studies, it is of the utmost importance 
to possess robust and advanced calculations based on concepts from traditional nuclear 
physics. Thereby, one of the major sources of uncertainty stem from the description of 
FSI mechanisms. In our eikonal model, we can either use optical potentials (ROMEA) or 
a Glauber framework (RMSGA). In some kinematic region of moderate hadron momenta 
both approaches can be used 28|. As they adopt very different underlying assumptions, we 
consider a comparison between the predictions of the two approaches as a profound test of 
the trustworthiness of either approach. We computed the transparency of the 4 He(7,£>7r~) 
reaction for kinematics at 6 c . m . = 70° and 90° with ejected proton momenta ranging from 
500 MeV/c to 1 GeV/c. As can be appreciated from Fig. [TU both descriptions yield a 
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similar shape, but the RMSGA calculations are consistently larger by about 5%. At higher 
nucleon momenta though, the difference is at the order of a few percent. 



a 
a, 

g 0.8- 



0.6- 



\t\ [GeV 2 ] 



FIG. 9: Contributions of the pion (dashed-dotted) and nucleon (dashed) to the total nuclear 
transparency (full) extracted from 4 He(7,p7r~) versus | t | at 6 c .m. = 70°. All calculations include 
CT. 




FIG. 10: The | t (-dependence of the relative increase of the nuclear transparency due to SRC 
and CT effects. We consider the 4 He(7,p7r~) reaction at # c .m. = 70° (left panel) and 90° (right 
panel) and kinematic conditions from Tabled The baseline result is the RMSGA calculation. The 
solid (dashed) curve includes the effect of CT (SRC). The dot-dashed line is the combined effect 
of CT+SRC. 



C. Pion electroproduction 

The E01-107 collaboration at Jefferson Lab has measured the nuclear transparency for 
the pion electroproduction process on H, 12 C, 27 A1, 64 Cu and 197 Au. Measurements were 
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FIG. 11: Comparison between the RMSGA (squares) and ROMEA (circles) description of the 
nucleon transparency of the 4 He(7,p7r~) reaction for kinematics at 9 c .m. = 70° (left panel) and 90° 
(right panel). Neither CT nor SRC effects were included in the calculations. 

done for the kinematics listed in Table [III In all the measurements the pion is detected in a 
relatively narrow cone about the momentum transfer. We have performed calculations for 
all target nuclei. The transparency is defined as 



(63) 



S a I dojY(u)) J A 3 pm d-Prn y J t\ -,/L T ./*» ,/<>\ / [.> MSG \ 



Ea / duY(u) J A3pm dp r , 



df! / dEi dE 7T (Kl 7V dQ,f 



lN / RPWIA 

The integration over uj takes into account the spread in energy of the virtual photon in 
the experiment and weighs each point with the reconstructed yield Y(oS) [42]. The quantity 
A 3 p m specifies the phase-space of the missing momentum and is determined by the condition 
| p m |< 300 MeV/c and the experimental cuts and detector acceptances. Accordingly, the 
final neutron is extremely slow and we have assumed that it possesses a transparency of 
one. A cut of 100 MeV was placed on the missing mass of the final state. The intranuclear 
attenuation effects on the ejected pion are again computed with the RMSGA model. We 
use a parametrization provided by the E01-107 collaboration for the free electroproduction 
cross section in Eq. (1461) 
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43]. 



Fig. [12] presents the results from our transparency calculations for the electroproduction 
reaction. The RMSGA calculations show a modest increase over the Q 2 range. This behavior 
finds a simple explanation in the p n dependence of the of Fig. [2J The results contained 
in Fig. [T2l cover a range in pion momenta given by 2.8 < p^ < 4.4 GeV. In this range, 
displays a soft decrease, which reflects itself in a soft increase of the nuclear transparency. 
The RMSGA+CT transparencies are again about 5% larger than the RMSGA ones. The 
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Q 2 E e 6 e E e , Pw 0, 



1.10 4021 27.76' 
2.15 5012 28.85' 
3.00 5012 37.77 
3.91 5767 40.38' 
4.69 5767 52.67 



1190 2793 10.58' 
1730 3187 13.44 
1430 3418 12.74 
1423 4077 11.53' 
1034 4412 9.09' 



TABLE II: Central values of Q 2 (GeV 2 ), incoming electron energy ^(MeV), electron scattering 
angle 6 e (degrees), scattered electron energy E e > (MeV), ejected pion momentum p w (MeV) and 
ejected pion angle (degrees) for the kinematics of the Jefferson Laboratory experiment E01-107. 
Angles are measured relative to the incoming electron beam. 

RMSGA+CT shows a strong Q 2 dependence with CT-related enhancements up to 20% at 
the highest energies. The evolution of the A-dependence of the transparency is shown in 
Fig. [131 One observes that the addition of CT to the calculation adds more curvature and 
that this increases with higher Q 2 . Finally, we compare our model calculations with the 



results from the semi-classical model of Ref. 



161 ]. The transparency is plotted as function 



of k = p-x — q. As in the photoproduction calculations [17| , our results again turn out to be 
higher by a few percent. 



IV. CONCLUSION 



We have outlined a relativistic framework to compute nuclear transparencies in exclu- 
sive A(7, Nir) and A(e, e'Nir) reactions. For the bound states, the model uses relativistic 
mean-field wave functions. At sufficiently high nucleon and pion energies, the intranu- 
clear attenuation on the ejected particles can be computed with a relativistic version of the 
Glauber model. At lower ejectile energies, the framework offers the flexibility to use opti- 
cal potentials. For nucleon momenta where both approaches can be applied, the Glauber 
and optical-potential based calculations predict nucleon transparencies in 4 He which follow 
similar trends. The differences in the magnitude of the transparency is smaller than 5% 
and shrinks with nucleon momentum. Our RMSGA predictions for the pion transparencies 
are in reasonable agreement with the semi-classical results of Larson, Miller and Strikman. 
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FIG. 12: (Color online) The Q 2 -dependence of the nuclear transparency for the A(e, e'7r + ) process in 
12 C, 27 Al, 63 Cu and 197 Au. The black and green curves are RMSGA and RMSGA+CT calculations 
respectively. The blue and red line are RMSGA+SRC and RMSGA+SRC+CT results. 

Both models predict similar trends, with the RMSGA predictions being systematically ~ 5% 
higher. This provides support that the baseline nuclear-physics transparencies can be com- 
puted in a rather model-independent fashion. Extension of our relativistic and quantum 
mechanical photoproduction calculations up to energies accessible in the JLab 12 GeV up- 
grade show an increase of the transparency up to 20% at the highest energies due to color 
transparency. Transparencies are also enhanced through the inclusion of SRC effects in the 
calculations. This yields an increase of about 5%, independent of the hard scale. Accordingly 
the SRC and CT mechanisms can be clearly separated. 



[1] A. S. Carroll et al., Phys. Rev. Lett. 61, 1698 (1988). 
[2] I. Mardor et al., Phys. Rev. Lett. 81, 5085 (1998). 

[3] A. Leksanov et al., Phys. Rev. Lett. 87, 212301 (2001), hep-ex/0104039. 



30 



1 2 3 4 56 10 20 30 100 200 

A 



FIG. 13: (Color online) A-dependence of the transparency for the A(e,e'ir + ) process at Q 2 = 
1.1 GeV 2 (black) and Q 2 = 4.69 GeV 2 (red). The solid curves denote RMSGA+SRC results. The 
dashed lines are RMSGA+CT+SRC calculations. 
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FIG. 14: Nuclear transparency results for C(e, e'7r + ) versus the z component of k = p n — q for 
kinematics corresponding to data points of the JLab experiment of Ref. [if]]. The circles are 
RMSGA+CT predictions, whereas the stars are from the semi-classical calculations of Ref. [it] ]. 
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